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SUMMARY 


The theoretical value of the integral of thermal conductivity is 
compared with the experimental values from shock-tube measurements. The 
particular case considered is the one-dimensional nonsteady flow of heat 
through air at constant pressure . This approach has Been previously 
described in NASA TR R-27. However, the correlation Between theory and 
experiment was uncertain Because of the large scatter in the experimental 
data. In this paper, an attempt is made to inprove the correlation By 
use of a more refined calculation of the integral of thermal conductivity, 
and By use of improved experimental techniques and instrumentation. As 
a result of these changes, a much closer correlation is shown Between 
the experimental and theoretical heat-flux: potentials . This indicates 
that the predicted values of the coefficient of thermal conductivity for 
high-temperature air may Be suitably accurate for many engineering needs, 
up to the limits of the test (^600° K) . 


INTRODUCTION 


Under normal conditions of temperature and pressure, the variation of 
the thermal conductivity with temperature for most gases is similar to 
that for the ideal gas, where the thermal conductivity is approximately 
proportional to the square root of the temperature . For these conditions, 
the only excited degrees of freedom are translational and rotational. 
However, the thermal conductivity deviates from the ideal-gas variation 
as the temperature increases to a point where the inactive degrees of 
freedom are excited and where chemical reactions are initiated (ref. l) . 
The chemical reactions cause the thermal conductivity to Become pressure- 
dependent as well. The mechanism of thermal conduction in a gas medium 
with a temperature gradient can Be explained as a transport of thermal 
energy. The excitation and chemical reactions enable the molecules to 
absorb thermal energy in high-temperature regions and to release the 
energy in low-temperature regions . 

In the case of air, the previously inactive modes which are excited 
are vibration of the molecules and transition of electrons to higher 
energy states . The important chemical reactions are the formation of 
nitric oxide, the dissociation of oxygen, the dissociation of nitrogen, 
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and finally, if a sufficiently high temperature is reached, the ionization 
of oxygen and nitrogen. Since each of these phases reaches its maximum 
at a different region of pressure and temperature, the variation of 
thermal conductivity is quite complex. 

Hansen (ref. 2) made calculations of coefficients of thermal conduc- 
tivity which included the aforementioned effects, with the exception of 
the formation of nitric oxide. In reference 3, the single effect of oxygen 
dissociation on the thermal conductivity of air was evaluated from 
measurements in a shock tube . The purposes of that experiment were to 
evaluate an average coefficient of thermal conductivity, and to check 
the theoretical calculations. The experimental results generally agreed 
with the theoretical estimates, though there was a rather large amount 
of scatter in the data. 

The purposes of this paper are: (l) to present a supplementary set 

of data which were obtained with improved experimental instrumentation 
and techniques; and (2) to present a more refined calculation of the 
integral of thermal conductivity for air in which the effects of the 
formation of nitric oxide have been included, and the major reactions are 
considered to occur simultaneously. 
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SYMBOLS 


number of mols of reactants (A) , also thermal diffusivity 

reactant of a chemical reaction 

number of mols of products, B 

product of a chemical reaction 

specific heat 

specific heat at constant pressure 

energy per mol at zero absolute temperature 

coefficient of thermal conductivity 

average coefficient of thermal conductivity 

coefficient of thermal conductivity for ideal gas with 
Sutherland correction (see eq. (24)) 

chemical equilibrium constant based on partial pressures 

nitrogen atom 
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n 2 

nitrogen molecule 

0 

oxygen atom 

0 2 

oxygen molecule 

p 

pressure 

p(Ai) 

partial pressure for component 

Q 

total partition function 

%> 

pQ, partition function for the standard state, 1 atm. pressure 

S (A i } 

partition function for component 

q. 

«T 

heat flow 

T 

temperature 

t 

time 

X 

distance 

X(Ai) 

mols of component A ± per mol of initially undissociated air, 
or mol concentration 

P 

density 

9 

r T 

/ k dT ; heat flux potential 
^0 

* 

q> 

P T 

/ k*dT, heat flux potential for ideal gas 
^0 

Subscripts 

• i 

ith component of gas mixture 

0 

condition at interface between hot gas and solid wall 

1 

condition in platinum film 
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THEORY 

Equilibrium Composition of Air 


The theoretical and experimental investigation of reference 3 on heat 
conduction in air covered a range of temperatures and pressures where the 
dissociation of oxygen and the formation of nitric oxide reach their 
maximums. However, when the coefficients of thermal conductivity were 
computed, only the effects of oxygen dissociation were treated because of 
the order of accuracy heing considered. The conductivity calculations of 
this paper have been refined by the inclusion of the effects of nitric 
oxide formation, nitrogen dissociation, and all inactive degrees of freedom. 
This is consistent with the improvement in the experimental procedures and 
techniques to be described in this paper. 

The general approach for the calculation of the equilibrium composition 
of air is: 

1. Establish from the available spectroscopic data, the partition 
functions of the various constituents of air (N 2 , 0 2 , NO, N, and 0) for 
translational, rotational, vibrational, and electronic excitations. 

2. Combine the partition functions to yield the equilibrium constant 
for each of the three reactions, oxygen dissociation (Os'* 20), nitrogen 
dissociation (n 2 -*• 2N) , and nitric oxide formation (N 2 + 0 2 -* 2N0) . 

3. Use the equilibrium constants to establish the set of simultaneous 
equations from which the equilibrium mol fractions are calculated. 

Consequently, the first step is to use the methods of statistical mechanics 
to formulate the partition functions for the various components in air. 

The detailed development of the use of partition functions in calculating 
the thermodynamic and transport properties of air is given in reference 2 . 

The partition function of a given molecule is defined as 
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which represents the sum of the exponential of every quantum level of 
energy for all states. The g± factor, or the degeneracy, accounts for 
the total number of states which have the same energy level, €j_. The 
magnitudes of the various energy levels, and the degeneracy of these 
levels have been accurately established from the spectroscopic data. 
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In a temperature range of 1000° to 6000° K and a pressure range of 
1 to 100 atmospheres, air may be assumed to be made up of N 2 > O 2 , NO, 
N, and 0. The resulting expressions for the natural logarithms of the 
partition functions of various species are 
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where the spectroscopic constants are obtained from reference 
The equilibrium constant, Kp, for each reaction, such as: 

SnpAp ZbpBp 


is then related to the partition function of the constituents of the 
reaction through the equation (ref. 5) 
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Zn Kp = _ prji + 
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where 


and 


Zn Qp = Zn Q + Zn p 


( 7 ) 
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AFi 0 — ZbpE 0 (Bp) — Z apE 0 (Ap) 


(9) 
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that is, the zero point energy of products less the zero point energy of 
reactants both referred to their standard states . Values of AE 0 for 
the various reactions considered here are found in Gilmore's report 
(ref. 6) . Substitution of these values into equation (7) results in the 
following expressions for the formation of nitric oxide, and the 
dissociation of oxygen and nitrogen 


Zn K r 


— 1 N2 + — O2 
2 2 


►NO 


10810 


59000 


+ Zn Qp(NO) - | Zn Qp(N 2 ) - | 0p(0 2 ) 

( 10 ) 

( 11 ) 


Zn Kp(0 2 -20) = — + 2 Zn Qp(0) - Zn 0p(0 2 ) 


Zn 


Kp(N 2 -2N) = - 11 ^ P °° + 2 Zn Q^N) - Zn Qp(N 2 ) 


( 12 ) 


According to reference 6, there is also nitric oxide dissociation; 

N0-*N + 0 (13) 

with a reaction energy between those for oxygen and nitrogen dissociation. 
As equation ( 13 ) does not introduce any element which the other reactions 
do not include, the number of mols for the different species (n 2 , 0 2 , NO, 
0, N) can be calculated with any three of the four reactions. Since at 
high temperatures oxygen and nitrogen dissociations are more important 
than nitric oxide dissociation, equations (10) , (ll) and (12) are chosen 
as the independent reactions for the subsequent analysis, and the reaction 
of equation ( 13 ) is treated as redundant. 
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The equilibrium constant, Kp, for each reaction, can also be expressed 
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through the definition of partial pressure 


p(a ± ) = 


X(A ± ) 


ZX(A ± ) + ZX(B i ) 


(15) 
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where X(A n - ) is the number of mols of species A± per mol of undissociated 
air. In the remainder of this discussion, X(Ai) will he referred to as 
the "mol concentration" of species Ap. The sum £X(Ai) + LX{B ± ) is the 
total mol concentration of both reactants and products, and p is the 
total pressure . 

In addition, the mol concentration of molecular nitrogen and oxygen can 
he expressed in terms of the number of mol concentrations of atomic 
nitrogen and oxygen, and nitric oxide by the following two balance equations 


X(0 2 ) = X^O*,) - | X(0) - | X(NO) 
X(N 2 ) = X 0 (N 2 ) - | X(NO) - \ X(N) 


(16) 

(17) 


where Xr,(0 2 ) and Xq(N 2 ) are reference quantities representing the constit- 
uents of one mol of air at one atmosphere and 293° K. The numerical values 

are: 


x q (n 2 ) = 0.78084 
x n (o 2 ) = 0.20946 


(18) 

(19) 


The remaining 0.0097 mol is composed of inert gases, such as A, He, etc. 
Substituting equations (15) , (l6) , (17) an d the constants of equations 
(18) and (19) into (l4) , one obtains the following three equations in 
terms of the equilibrium constants: 

[4 - Kp 2 (N0)] X 2 (N0) - Kp 2 (N0) [X(N0)X(N) + X(N)X(0) + X(N0)X(0) 


_ 1.9806 X(H0) - 0.41892 X(N) - 1.56168 X(0) + 0.65422] = 0 

( 20 ) 


4 + 


Kn(0) 


x 2 (o) + 


K p (0) 


X(N0)X(N) + X(N)X(0) + X(N0)X(0) 


I.9806 X(H0) - 0.41892 X(N) + 1.56168 X(0) - 0.82971 


= 0 


( 21 ) 
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Kp(N) 

p 


x 2 (n) 


K p (N) 

P 


X(NO)x(N) + X(N)X(0) + X(HD)X(0) 


+ I.9806 X(NO) + 0. 41892 X(N) - 1.56168 X(0) - 3.09306 


= 0 


( 22 ) 


The equilibrium mol concentrations of the various species were calculated 
by solving equations (20) , (2l) , and (22) simultaneously. The computation 
was done on an IBM 70^ computing machine and the results are illustrated 
in figure 1 as constant-pressure functions of temperature . Values for 
NO found in this manner were within 1 percent of the results of the 
detailed calculation performed by Gilmore (ref. 6) . The mol concentrations 
of NO, 0, and N are compared in figure 2. 


Thermal Conductivity 


The thermal conductivity of high-temperature air based on the number 
of mols of each species is now calculated by taking into account the 
effects of vibrational and electronic excitation as well as chemical 
reactions . It is found in the present analysis that the maximum value of 
mols of NO in air at pressures up to 100 atmospheres is only 11 percent 
(fig- 1( a) ) , and it has been observed that the NO molecules have vibra- 
tional characteristics similar to those of O 2 and N 2 (ref. 4) . Hence, 
the vibrational and electronic contributions to thermal conductivity in 
air may be found very closely if it is assumed that there are N 2 , O 2 , 

N and 0 only, and that there is no NO species . Such a procedure is 
followed by taking directly those values in reference 2 for the vibrational 
and electronic part of thermal conductivity. The other part of thermal 
conductivity due to chemical reactions may be obtained by the method of 
Brokaw and Butler (ref. 7 ) for simultaneous reactions. 1 The rigorous 
calculation is rather lengthy and complicated. Fortunately, in the case 
of simultaneous reactions, such as NO formation and O 2 dissociation, 
the over-all thermal conductivity can be approximated closely by the sum 
of thermal conductivities due to each individual reaction (see appendix A) . 


1 The collision cross section for diffusion is needed for the Brokaw 
and Butler formula. As previously pointed out by many authors, the 

collision cross section for diffusion is usually different from that for 
heat conduction and viscosity. In the case of the NO formation reaction, 
no accurately determined data exist for the collision cross section of 
NO with O 2 and N 2 - Hence, it is assumed that the collision cross 
section for diffusion is 0.778 of that for conduction and viscosity, that 
is, a molecule with a l/5-power repulsive force (ref. 8, p. 173 ) • 
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The expression for the total coefficient of thermal conductivity can 
then he expressed as: 

k total = k * + k vE + k N0 + ^ + k N ( 2 3) 


+ 

) 

L 


where k is defined as the coefficient of thermal conductivity of air 
which has no vibrational degree of freedom, dissociation, or ionization. 
Mathematically, it is expressed in reference 2 as: 

k* = 1.994X10- 3 JouleB (24) 

1 + cm-sec-°K 

T 

The term ky-g is defined as the coefficient of thermal conductivity due 
to electronic and vibrational excitation of all molecules, kQ, the contri- 
bution due to oxygen dissociation, kjj, the contribution due to nitrogen 
dissociation, and k^o* the contribution due to nitric oxide formation. 

The component thermal conductivities are shown in figure 3 in terms of 
the ratio of the coefficient k to the reference quantity k*. 


Comparisons of the partial thermal conductivities for the various 
degrees of freedom are shown in figure The effect of molecular 
vibration becomes appreciable at temperatures above 500° K, that of NO 
formation above 1500° K, and that of O 2 dissociation above 2000° K. 

It may also be observed that the slope of the curve is greater if the 
activation energy involved is greater. Thus, the curve for vibrational 
and electronic excitation has the smallest slope whereas the curve for 
0 2 dissociation has the largest slope. The thermal conductivity for the 
various degrees of freedom becomes pressure dependent toward higher 
temperatures, as seen from the separation of the curves for the various 
pressures. The separation of the curves for the combination of vibrational 
and electronic degrees of freedom is due to the pressure dependence of the 
formation of atoms which become electronically excited at higher temper- 
atures . The separation of the curves for NO formation and O 2 disso- 
ciation are reflections of changes in the equilibrium mol concentrations 
of species involved in the reactions. In general, the pressure dependence 
of the thermal conductivity of air can be attributed mainly to oxygen 
dissociation. For the range of temperatures investigated the total 
coefficient of thermal conductivity is shown in figure 3(e) • 

Heat Conduction Between Two Semi-Infinite Media 

The computed value of thermal conductivity can now be related to the 
specific case of time-dependent one-dimensional heat conduction through 
a semi-infinite gas medium at constant pressure with uniform initial 
conditions. The partial differential equation governing this process is 
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dx 
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This equation can he solved hy specifying the quantity 


(25) 


T 

? - [ *total dT - & (26) 


as the dependent variable instead of T itself (ref. 3)* This quantity 
will he called the heat-flux potential. Substitution of equation (26) 
into equation (25) results in the following partial differential equation 
with cp as the dependent variable: 

69 _ a ^JP = o (27) 

dt 6x 2 


where the thermal diffusivity, a = k/CpP. 

The theoretical components of cp are obtained hy integrating the 
theoretical values of each component of k in equation (23) according 
to equation (2 6) . The total heat-flux potential is expressed in dimen- 
sionless form hy dividing each component hy cp* • This quantity is defined 
as the heat-flux potential for ideal air and is given hy: 


T k* dT (28) 

where k* is given hy equation (2^) . Numerical values for this integral 
are shown in figure 5 • 

The total heat-flux potential, expressed in dimensionless form, cp/ cp* , 
is then equal to 
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vCpVVE 


-cp /NO 




(29) 


where the 1.0 denotes the idealized air value, (cp/cp )yg the contributions 

from vibration and electronic excitations, ( cp/cp* ) N0 the contribution of 
nitric oxide formation, (cp/cp*) 0 the contribution of oxygen dissociation, 
and ( cp/cp*) N the contribution of nitrogen dissociation. Variations of 
these components with temperature are shown in figure 6. 

The heat-flux potential as a function of temperature can also be 
obtained by combining the solution of the partial differential equation (27) 


H o *=- > 


<£j O i — I 
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with experimental measurements of certain "boundary conditions. The 
equation can be solved if two boundary conditions are available: (l) the 

interface temperature at the wall which is evaluated by the use of a thin 
film gage (a detailed description of gage construction and performance is 
given in the experiment section) , and (2) the initial state of the high- 
temperature air in the region behind the reflected shock, which is a known 
function of the incident shock-wave speed. The measurement of these two 
quantities permits the calculation of the gradient of thermal conductivity 
which exists in the gas at the interface, where the temperature rise is 
relatively small. Once the value of the gradient at the interface is 
specified, one can integrate the differential equation (27) to obtain the 
value of the integral of thermal conductivity for the high-temperature air 
in the region behind the reflected shock. The qualitative form of the 
solution for the limiting value of cp as a function of (dcp/ox) 0 is shown 
in sketch (a) . The numerical value of cp is determined from these graphi- 
cal solutions using the measured value 
for (dcf/dx) 0 * Details of this approach 
can be found in reference 3. These 
solutions will be compared with the 
theory in a subsequent section of this 
report . 


EXPERIMENT 

Only a brief description of the 
shock-tube apparatus is presented in 
this paper since a detailed descrip- 
tion is available in reference 3* 
This section will be more concerned 
with improvements in experimental 
procedures, and improvements in the 
design and manufacture of the film 
gages used to measure temperature 
at the wall of the shock tube . 


* 



Sketch (a) 


The shock tube used to generate high-temperature air for the experi- 
ments consists of a reservoir chamber initially filled with high-pressure 
helium and a test tube initially filled with air at relatively low pres- 
sures (usually less than one atmosphere) . The two sections are separated 
by a diaphragm. Upon rupture of the diaphragm, a strong shock forms near 
the diaphragm, travels down the tube at high speed, accelerating and com- 
pressing the air after it, and then on reaching the end of the tube is 
reflected back toward the reservoir. The actual test phase begins at the 
instant the incident shock wave is reflected from the end of the shock 
tube . At this instant the air sample has been brought to rest by the 
reflection so that heat diffusion in air which is in a state of high tem- 
perature and pressure can be studied in the absence of strong heat con- 
vection effects . Relaxation effects do not seem to be large since the 
relaxation time is of the order of 10 microseconds as compared with the 
minimum testing time of 200 microseconds . Temperature detection is accom- 
plished by a thin film resistance gage mounted flush with the wall at the 
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end cap of the shock tube . A schematic diagram for the test section is 
shown in sketch (b) . During the test period, the reflected shock wave 

Pressure 



travels faster than the cooling influence of the heat -diffusion process . 
Moreover, the heat penetration into the pyrex backing of the film gage 
is so small that the l/4-inch thick pyrex is practically an infinite heat 
sink throughout the test. Hence, there occurs the one-dimensional, 
nonsteady heat conduction described in the Theory section. 

The partial differential equation (eq. (27)) governing this process 
requires two boundary conditions for its solution. Both are found from 
the shock-tube experiment. The first boundary condition required is the 
initial state of the high -temperature air. The shock-wave speed is 
determined by measuring the time for the shock to traverse two pressure 
pickups spaced at a known distance . The state of the air in the test 
region is then obtained from reference 9 which presents the state variables 
behind a reflected shock as a function of the shock speed. The second 
boundary condition, the interface temperature at the wall, is measured by 
the thin film- resistance gage. (A discussion of this thin film gage is 
presented in appendix B.) This measurement is the most critical one of 
the entire experiment, as it is believed that the scatter of the data in 
reference 3 can be attributed mainly to the erratic behavior of the film 
gages used in that experiment. The improved gage used for this experiment 
is of the same type, a resistance thermometer with a microsecond response 
time. The measuring circuit is a Wheatstone bridge with the film gage as 
one of its arms. The signal is recorded on an oscilloscope. 

The construction of the improved gage is somewhat different from that 
described in reference 3. The gage which yielded the most satisfactory 
results for the present series of tests is shown in figure 7(a) in four 
stages of manufacture. These stages are: 

1. A rectangular prism with grooved sides is formed from pyrex to 
serve as a base for the gage. 

2. The grooves are coated with three layers of Hanovia Platinum- 
Bright solution no. 5* The gage is baked after each coating at 1100° F 
for 20 minutes, followed by a cooling period of 8 hours or more. 
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3 . The film gage is formed by painting a platinum strip l/l 6 inch 
wide across the face of the prism from one groove to the other to form 
a continuous path. The gage is baked in the same manner as in step 2. 

4. Lead wires are soft soldered to the coated grooves. Precautions 
are taken against overheating during the soldering process so as not to 
destroy the resistance characteristics of the film. The purity of the 
soft solder is also checked as it was found that impurities increased the 
terminal resistance and decreased the gage sensitivity. 

Figure 7 (b) shows the completed gage mounted in the holder with a 
resin cement. The final step is to coat the face of the gage with a 
1000 A layer of electrically inert silicone monoxide. The purpose of 
this coating is to protect the film from free electron or ion interference 
from the partially ionized air. Initial experiences with uncoated film 
gages showed spurious, fast pulse signals of large amplitude at higher 
temperatures where ionization occurs. Similar experiences were reported 
in reference 10 and were attributed to charge pickup from the ionized gas. 
It is assumed that the coated gages, which produce outputs free of such 
spurious pulses, indicate the correct temperature except for the time 
lag effect produced by the silicone monoxide coating. 

The experimental procedure was also changed from that used in the 
experiments of reference 3 . The following steps also tended to eliminate 
some of the scatter in the data: 

1. The original film gages were calibrated in a water bath 
calorimeter (see ref. 3 ). However, this type of calibration is very time 
consuming and does not simulate the fast step impulse of temperature 
obtained under actual test conditions. A capacitor discharge calibration 
was tried, similar to that reported in reference 11 , but this also did 
not simulate the condition that heat flux varies as the square root of 
time. For a true simulation, the shock tube itself was used for calibra- 
tion of the film gage at gas temperatures below 1000° K (i.e., low shock 
strength) . At these test conditions the interface temperature jumped 
approximately 10° K above atmospheric temperature after the reflection of 
the shock wave. At these temperatures where the thermal conductivity has 
been established, an equivalent constant for temperature resistivity was 
found by equating the experimental conductivity to the known values. 

This constant was then used to analyze the data for higher temperature 
tests where the interface temperature rise varied from 20° K to 60° K. 

2. Prior to the test, both the reservoir chamber and test tube were 
evacuated to eliminate any contamination. In addition, the test tube 
underwent a purge cycle three times. 
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3. During testing, a high-frequency vibration produced by the rupture 
of the diaphragm, was evidenced by a fluctuating signal from the film 
gage . Further examination revealed that the vibration was being trans- 
mitted to spring-type connections at the open ends of the lead wire . A 
solder connection at this location eliminated the spurious signal. 

Some typical temperature records registered by these gages are shown 
in figure 8- The appropriate test conditions for each record are 
indicated on the figure. 


DISCUSSION OF RESULTS 


The theoretical values of cp/cp* are obtained from figure 6 for the 
experimentally determined temperatures and pressures. The variation of 
pressure with increasing temperature (i.e., increasing shock strength) 
for this series of experiments can also be seen in figure 6 (b) . The 
over-all potential is obtained theoretically by evaluating the components 
cp/cp* for the various modes of excitation and the various reactions, then 
algebraically adding the components. At temperatures ranging from 600° K 
to 2000° K, the increase in the over-all heat-flux potential above the 
ideal-gas value of 1.0, can be attributed to the vibrational and electronic 
modes of excitation. The formation of nitric oxide starts at 1000 K, 
although its effect on cp/cp* does not become noticeable until a tempera- 
ture of 2000° K is reached. At this point, oxygen dissociation is ini- 
tiated, and by the time a temperature of 3000° K is reached, the cp/cp 
contribution due to oxygen dissociation is of the same magnitude as those 
for the vibrational and electronic modes and for nitric oxide formation. 
Above 4000° K the theoretical value of the dimensionless heat-flux poten- 
tial levels out slightly and then increases again. This variation 
indicates that the nitric oxide formation and oxygen dissociation are 
nearly complete, and that nitrogen dissociation is beginning. 

The theoretical value of cp/cp* is compared in figure 9 with the 
values of cp/cp* obtained by using experimental boundary values to solve 
the one-dimensional heat-flow equation (eq. (23)). The experimental and 
theoretical cp/cp* agree very well in a temperature range from 1000° to 
4600° K. Since the coefficient of thermal conductivity, k/k*, can be 
derived from the slope of the temperature versus heat-flux curve (fig. 9) > 
the comparison of experimental and theoretical heat-flux potentials is 
then a direct indication of the accuracy with which the coefficient of 
thermal conductivity has been predicted. However, the uncertainty in 
the determination of a derivative from any experimental data is relatively 
large. In this light, the close correlation between experimental and 
theoretical estimates of the heat flux can be interpreted as a necessary 
but not sufficient criterion that the coefficient of thermal conductivity 
has been accurately predicted. The close correlation does indicate that 
the predicted coefficients of thermal conductivity may have suitable 
accuracy for many engineering needs. 
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CONCLUDING REMARKS 


Computed values of the integral of the coefficient of thermal conduc- 
tivity for high- temperature air were checked hy Hansen, et al., in NASA 
TR R-27, hy comparing the theoretical values of the integral of thermal 
conductivity (heat-flux potential) with the experimental values of the 
integral obtained from shock- tube measurements . However, the correlation 
between theory and experiment was uncertain, mainly because of the large 
scatter in the experimental data. A much closer correlation has been 
obtained in this paper by using a more refined calculation of thermal 
conductivity, and by using improved experimental techniques. 

The theoretical integral of thermal conductivity is obtained by 
algebraically adding the integrals for the excitations of the various 
degrees of freedom for the components of air, and to the chemical 
reactions between these components. This calculation must include the 
excitation of the vibrational and electronic modes, and the formation of 
nitric oxide, as well as the dissociation of oxygen considered in 
reference 3- The inclusion of these additional effects is necessary 
because they determine almost entirely the deviation of thermal conductivity 
from the perfect gas value for temperatures up to 3000° K. 

The most critical measurement of the entire experiment is the 
measurement of the interface temperature at the wall of the shock tube 
by a thin film resistance gage. Consistent and accurate temperature data 
are obtained if the following steps are taken: (l) positive binding between 
the painted platinum film strip and the pyrex backing; (2) protection of 
the film and terminal connections from extraneous electrical pickup and 
damage from the high-temperature air, and (3) calibration of the resistance 
gage during actual shock-tube experiments at low temperatures where the 
variation of the coefficient of thermal conductivity with temperature has 
been definitely established. 

As a result of these factors, a much closer correlation is shown 
between the experimental and theoretical integrals of thermal conductivity. 
This close correlation is a necessary but not sufficient criterion that 
the coefficient of thermal conductivity has been accurately predicted. 

The close correlation does indicate that the predicted coefficients of 
thermal conductivity may have suitable accuracy for many engineering 
needs, up to the limits of the test (4600° K) . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Oct. 31, i 960 
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APPENDIX A 

♦ 

INDEPENDENCE OF THERMAL CONDUCTIVITY 


In a pure dissociating gas, such as 0 2 or N 2 , the thermal conduction 
at high temperatures can he ascribed to two simultaneous processes: (l) 

conduction by molecular and atomic collisions, (2) conduction in the form 
of transported chemical enthalpy of dissociated molecules due to concen- 
tration gradients (i.e., chemical diffusion). For a gas mixture such as 
air, where both of these species are dissociating and NO is forming 
and dissociating, the basic mechanism for thermal conduction is the same. 
Collision occurs between particles of all the components in the mixture 
(i.e., 0 2 , NO, N, etc.). However, the transport rate of chemical 
enthalpy of dissociated molecules for each reaction is affected by the 
presence of the other reactions, though the other reactions may be 
quiescent at a given state condition. This effect is included by Butler 
and Brokaw (ref. 7) in the following expression for the thermal conductivity 
due to chemical diffusion. 

0 AH3 . ah 2 

AH3. An A 12 

AH 2 A 2 1 A 22 

An Ai 2 

A 2 1 A 22 

for the case where two reactions are simultaneously occurring. The AH f s 
are the heats of reaction for each reaction alone, A^ is a measure of 
the reaction rate of the ith reaction alone, and Aj_j is a measure of 
the change in the reaction rate of the ith reactions due to the presence 
of the jth reaction. The A T s are expressed in reference 7 as functions 
of the number of mols and the collision cross sections for the various 
combinations of reacting particles . 

Rewriting the expression for k r , and setting 

Ai 2 
A22 

^ AHi 2 A 22 AH 2 2 Ah 

+ A 2 i) — — — 

|a| IaI J 

A12A21 « AhA 22 , 


I A I = 


An 

A21 


one has, 


kr = 


RT" 


AHiAH 2 

IaI 


(A 


12 


Now, if A i 2 = A 2 i An or A 22 , and 



■t 


f-* Q -P- > 


2T. 


17 


then. 


k T = — + 1 ^ 2 ‘ 


RT 2 An 


RT 2 A 22 


The first tern on the right side of the above equation is the thermal 
conductivity due to the first reaction alone as indicated by AHx, and, 
similarly, the second term is due to the second reaction with AH 2 only. 
Hence, for the approximations stated above, the over-all thermal conduc- 
tivity, k r , is equal to the sum of thermal conductivities that result 
from first and second reactions, respectively. 

The magnitudes of i\ij T s for the simultaneous reactions will now be 
examined. The atomic collision cross sections for diffusion were 
obtained from reference 2, and the molecular collision cross sections for 
diffusion were obtained by correcting the collision cross sections for 
conduction in reference 2 by a factor of 0.778 as suggested in reference 
8 for a molecule with a l/5-power repulsive force. In the case of three 
simultaneous reactions, such as NO formation, 0 2 dissociation, and 
N 2 dissociation, the N 2 dissociation has only a small effect in the 
region of temperature considered in this report; and between NO formation 
and 0 2 dissociation, the magnitudes of the Aj[j terms are: 


*12 


= A 21 = 2.681X10 5 


An = 29.91X10 
A 22 = 22 . 44x10 5 


cm-sec 

mol 


cm-sec 


mol 


cm-sec 




XL2 „ 1 


tLl 


10 


mol 


A 12 1 

A 22 10 


for T = 3348' k 
p = 10.74 atm 


A 12 = A 2i = 2.389x10 s 

cm-sec 

mol 

An = 66.51X10 5 

cm-sec 

mol 

a 22 = 1.029x10 s 

cm-sec 

mol 

A 12 ~ 1 A 12 p 

An 30 A 22 500 





for T = 2647° K 

p = 19*69 atm 


l8 


These calculations show that it is reasonable to assume that the 
over-all thermal conductivity due to WO formation, 0 2 dissociation, 
and W 2 dissociation is equal to the sum of the contributions from 
each reaction. 

In a private communication, Dr. Brokaw stated that where the 
off-diagonals (i.e., A 12 ) are as much as 20 to 50 percent of the diagonal 
terms (i.e., A 1X , A 22 ), the error in thermal conductivity introduced by 
neglecting the off-diagonals is only 5 percent at worst. 




APPENDIX B 


THEORY OF THIN FILM THERMOMETRY 


The theory of thermometry for a thin film gage may he outlined as 
follows: Consider a one-dimensional nonsteady heat conduction (sketch c) . 

Heat is being transferred from the hot 
q Hot gas gas through a layer of metal to glass. 

It is physically true that the tempera 
ture and heat flow on the gas side of 

l 1 1 ( metal are different from those on the 

glass side. The difference is caused 

777777/77777777777777777mT77 by the fact that the metal needs heat 

to adjust its thermal level with the 
gas and the glass on two sides and 
that it takes a finite time for the 
heat-conduction process to be complete 
However , both temperature and the heat-flow differences can be reduced to 
a negligibly small amount for a thin gage so that the heat-conduction 
process will be maintained practically the same as if there were no metal 
film in between. 


Pyrex backing 

Sketch (c) 


It is further reasoned from reference 3 that for a gas instantly 
heated in a shock tube, the interface temperature between the metal and 
gas will remain constant and the heat flow at that interface will be 
inversely proportional to the square root of time. Mersman, Berggren, 
and Boetler arrived at essentially the same conclusion in dealing with 
heat conduction between two semi-infinite solids (ref. 12). If only the 
heat conduction inside the metal is considered and a perfect contact is 
assumed, the temperature and heat flow at any point inside the metal may 
be obtained from the equations in reference 12. Using temperature, T Q , 
and heat flow at the interface, Qq, as boundary conditions, one obtains 
the following expressions: 


T(x,t) 

T 

^O 

Q(x,t) 



X 2 


e 4ait 


where 


- T, 
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where Tj_ is initial temperature of metal and a^, is the thermal 
diffusivity of the metal, which is taken to he a constant . For 
x = 2000 A, t = 1 psec, and a platinum film, 


T(x,t) 

T 

- L o 


l - 0.0027 = 0.9973 


■ - X — = 1 - 2.4x10 4 = O.99986 
^o 

The results are applicable for the case where there is no glass hacking, 
hut a semi-infinite extension of metal. Since glass has smaller thermal 
conductivity, the deviations for temperature and heat flow will certainly 
he much less than indicated above. Hence, one can he assured that in 
using the thin film gage of only 2000 A thickness, one obtains essentially 
the surface temperature and heat flow for the hacking material. 


<-=r o h 


<: -4- o h 
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Mol concentration of nitric oxide, X(NO) 
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(a) Nitric oxide. 

Figure 1.- Equilibrium mol concentration of air components. 
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Figure 1.- Continued. 
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Mol concentration of atomic nitrogen, X(N) 
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(c) Atomic nitrogen. 


Figure 1.- Concluded. 
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Figure 2.- Component mols of air in equilibrium. 
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Partial coefficient of thermal conductivity. 
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(a) Nitric oxide formation. 

Figure 3 *- Coefficients of thermal conductivity. 
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(To) Oxygen dissociation 0.1 to 0-9 atm. pressure. 
Figure 3*- Continued. 
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(c) Oxygen dissociation, 1 to 100 atm. pressure. 
Figure 3 • - Continued. 
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(d) Nitrogen dissociation. 


Figure 3 -~ Continued. 
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(e) Total coefficient of thermal conductivity • 
Figure 3 • - Concluded . 
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Figure 4.- Comparison of the partial coefficients of thermal conductivity. 
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Figure 5.- Integral of reference thermal conductivity. 
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(a) Partial heat flux potentials. 
Figure 6.- Heat flux potentials. 
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(t>) Total heat flux potentials. 


Figure 6.- Concluded. 
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(h) Gage mounted in holder. 

Figure J.- Photographs of thin film gage. 
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(a.) M]_ — 5*88, P3 = 12.96 atm, 
T 3 = 3294° K 
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(l) M x = 7-02, p 3 = 7-72 atm, 
T 3 = 3956° K 
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(a) Mi = 8.525, P 3 = 1-76 atm, 
T 3 = 4630° K 
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(e) Mi = 2.079, P3 = 23.I atm, 
T 3 = 686.5° K 


Figure 8.- Records of temperature at the surface of a pyrex slab. 
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Figure 9.- Comparison of experimental and theoretical integrals of thermal conductivity. 
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